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We study forward-backward charge fluctuations to probe the correlations among produced par- 
ticles in ultra relativistic heavy ion collisions. We develop a model that describes the forward- 
backward dynamical fluctuations and apply it to interpret the recent PHOBOS data. Within the 
present model, the dynamical fluctuations are related to the particle production mechanism via clus- 
ter decay and to long range correlations between the forward and backward rapidity hemispheres. 
We argue that with a tight centrality cut, PHOBOS may see a strong decrease of the dynamical 
fluctuations. Within the present model, this deterioration of the correlation among the produced 
hadrons can be interpreted as a sign for the production of a hot, dense and interacting medium. 
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I. INTRODUCTION 

One of the main goals of the heavy ion program is to 
understand the nature of the hadron production mech- 
anism (e.g. parton coalescence, string fragmentation or 
cluster decay) . Recent RHIC data on jet quenching and 
elliptic flow («2) can be interpreted as an evidence of 
forming a quark gluon plasma (QGP) during the colli- 
sion of two heavy gold nuclei at center of mass energy 
^/JWW = 200 GeV Q. However, the detailed mecha- 
nism how the deconfined partonic matter transforms to 
hadrons is still unknown. 

Using correlations and fluctuations to probe the nature 
of the created QCD matter has been pro pos ed b y many 
authors 
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[H E fia 12011211 l2aTl23| . Unfortunately, the experi- 
mental exploration of most of the suggested fluctuation 
signals did not yield positive evidence for the formation 
of a QGP. Nevertheless, this must not be seen as a con- 
tradiction to the formation of a partonic state at RHIC 
or SPS energies, because the proposed signatures may be 
weakened due to many effects (e.g. the phase transition 
is not a strict first order transition - which is the case 
for RHIC energies). Other possibilities are that rescat- 
tering and/or re-thermalization in the late QGP or in 
the subsequent hadronic stage may blur the fluctuation 
signals. Furthermore, the fluctuation signals might be 
modified and cloaked by the hadronization process it- 
self. However, we want to point out that recent studies 
of rescattering and/or diffusion in the hadronic environ- 
ment, seem to indicate that the diffusion in the hadronic 
stage might be small enough to allow a survival of the 
initial correlations |24j. 

In this letter we explore the forward-backward charge 
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fluctuations j2j| |2|| 153, to probe the particle produc- 
tion mechanism via cluster decay. Long ago, the number 
k of charged particles produced per cluster decay was 
estimated in pp and pp interactions by the UA5 exper- 
iment 29]. It was found that the effective multiplicity 
per cluster is k ps 2. Recently, the PHOBOS experi- 
ment performed a similar analysis for Au+Au reactions 
at ^/snn — 200 GeV and reports that k = 2.7 for pe- 
ripheral collisions and k = 2.2 for central collisions. Note 
that all measured cluster multiplicities are larger than 
expected for a hadron resonance gas (&hg = 1-5) Q, in- 
dicating that the measured charge correlations can not be 
described by simple statistical models based on hadronic 
degrees of freedom. 

Our goal is to explore how the cluster structure changes 
if deconfined matter is formed. The meaning of the 
term cluster is rather general, for example it could be a 
hadronic resonance, a partonic string or a QGP droplet. 
Either the survival or the destruction of the clusters pro- 
vides indirect information about the properties of the sur- 
rounding medium. For example what is the percentage 
of particles produced by cluster decays, are the hadronic 
expectations compatible with the data or not and if a hot 
and dense (partonic?) medium is formed or not. 

This letter is organized as follows: In Section [n] we 
develop a model that describes the multiplicity fluc- 
tuations 0% as measured by the PHOBOS experiment 
Eel I27L l28j . This measure is sensitive to charged 
particle production via cluster decay. In Section IIIII we 
apply our model to interpret the PHOBOS data and to 
estimate the fraction / of particles produced via clus- 
ter decay. In Section IIVI we discuss the results of our 
analysis of the measured data [2{J [2(| [2?], for central 
Au+Au collisions. Next, we address how the formation 
of a dense medium affects the dynamical fluctuations be- 
cause it might influence both the cluster multiplicity and 
the short range rapidity correlation length. Finally, we 
comment on the behavior of the dynamical fluctuations 



2 



in central collisions. In Section[V]we summarize the main 
results of this letter. 



contribution to the total fluctuations can be calculated 
from 



-7J+A77/2 jj+Atj/2 



II. THE MODEL 

In this section we develop a model to describe forward- 
backward charge correlations. We define two symmetric 
rapidity regions at ±77 with equal width A77. The num- 
ber of charged particles in the forward rapidity interval 
77 ± A77/2 is Nf while the corresponding number in the 
backward hemisphere —77 ± A77/2 is given by Nb- We 
define the correlation variable C such that 



C = 



N F -N 



B 



(1) 



The event-by-event fluctuations (variance) of C are 
decomposed into two parts, weighted by the respective 
fractions / and (1 — /) of particles from the different 
sources: 
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The first term in Eq. J2J cr| R denotes the short range cor- 
relation due to cluster decays. This term is a direct mea- 
sure of the particle production via cluster decay. As dis- 
cussed above, the specific nature of the clusters is rather 
general. 

°SR depends on the observability of the decay daugh- 
ters in the rapidity width A77, so that when A77 is larger 
than the size of the cluster one is able to observe all the 
decay products of this cluster. However, if A77 is small 
one may "loose" the decay products of such a cluster. To 
model this rapidity window effect we propose the follow- 
ing form 



<7g R = k [1 - exp(-A77/A sho rt)] 



(3) 



where A s hort is the rapidity correlation length between 
particles produced as a result of the cluster decay and k 
is the aforementioned multiplicity per cluster. 

The second term u£ R in equation (01 includes the fluc- 
tuations due to long range correlations between particles 
in both hemispheres and the background particles that 
are statistically independent. To describe the long range 
correlation, it is convenient to introduce a two body dis- 
tribution function P2 (?7i , ^2) ■ Such that 



= 1-A / d m 



dr}2 exp 



-ri-Ari/2 rj-Arj/2 



ZA long 



(5) 

Next, we decompose the total multiplicity N in each ra- 
pidity window A77 in two parts such that N — iVg R -|- ./Vl R 
where iVgR is the charged particle abundance produced 
from the decay of all clusters inside A77 and -/Vlr is the 
number of particles migrated from outside into the ra- 
pidity window under investigation. 

In each rapidity bin A77, the fraction of particles pro- 
duced by cluster decay / is given by / = Nsr/N. With 
the weighting of <r§ R and cr LR by / and (I — /), and 
the abbreviation £ = (1 — f)A one rewrites Equation © 
using Eqs. © and JSJ) as 
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where -F(Ai ong , A77, 77) is the result of the double integral 
in Equation (J5J. Equation (jfJJl includes different physical 
situations, e.g. in case of vanishing short range correla- 
tion / = 0, and Uq is determined solely by long range 
correlations. In case of / = 1, only short range corre- 
lations (i.e. cluster decays) contribute to the dynamical 
fluctuations. Therefore, if measurements are done such 
that the two bins centered around ±77 are far from each 
other, the contribution from long range correlations be- 
comes negligible. Only in this case, an increase of the 
bin size A77 to A77 3> A sn ort leads to the commonly used 
approximation ~ k from which the true cluster mul- 
tiplicity can be deduced. Another limiting value obeyed 
by Eq. © is that for a very large 77 the dynamical fluc- 
tuations tip do approach the poisson value er^ = I. This 
is due to the vanishing of the long range correlations and 
the low multiplicity which drives the fraction of corre- 
lated particles to be zero. 

Before we analyze the PHOBOS data let us shortly 
summarize the assumptions used in the above deriva- 
tions: (I) There is only one type of clusters. (II) All 
clusters decay to the same number of particles k, i.e. 
fluctuations of k itself are not included in the present 
model. 
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where A is the amplitude of the correlation function, 
giving the strength of the two body distribution when 
Vi = V2- Generally, A is a product of the true ampli- 
tude and a function that might depend on the motion of 
the center of mass rapidity of the pair ?7 C = (771 + 772 ) / 2 . 
For the present study, we assume that A is independent 
of 77 and A77. Under these assumptions, the long range 



III. RESULTS 



In this section we study the dynamical fluctuations 
for Au+Au collisions at ^snn = 200 GeV. We extract 
the values off and £ from PHOBOS data reported in 
[iJ|2l|2l|2||. Tabled summarizes our findings on the 
parameters of the present model. With these parameters, 
we predict a^, as a function of 77 and At;. Figure ^ de- 
picts both PHOBOS data for 0% - 20% central Au+Au 
collisions (symbols) and our model results (line) for the 



3 



data 


/ 


£ 


^short 


Along 


Au+Au 0% - 20% 


0.88 


1.8 


0.4 


0.7 


Au+Au 40% - 60% 


0.99 


2.5 


0.6 


0.9 



TABLE I: Parameters of our model as estimated by analyzing 
0%-20% central and 40%-60% peripheral Au+Au PHOBOS 
data at y'snn = 200 GeV using equation JSJ. 




FIG. 1: (Tc versus rapidity rj for 0% - 20% central Au+Au 
at A 77 = 2. Black squares are PHOBOS data 0, M El E3 . 
The line is the model calculation using equation @. 



dynamical fluctuations Oq as a function of r\ while the 
width of the observation window Ar) is kept as Arj = 0.5. 

In Fig. |21 we show cfq versus Arj, while the center of 
the observation window is fixed at rj — 2. We emphasize 
that the same parameter set (see TableHJ is used for both 
data sets in Figs. ^andEl 

Before we discuss the interpretation of our model re- 
sults, let us turn to semi-peripheral collisions. Figs.|3|and 
H depict the PHOBOS results on 40% - 60% peripheral 
collisions. As for central collisions, the model allows to 
describe both data samples with the same parameter set 
simultaneously (cf. Table P). The dynamical fluctuations 
o-q measured at rj = 2 for 40% — 60% peripheral Au+Au 
collisions and a very wide rapidity windows Arj = 2 yield 
an effective cluster multiplicity of k sa 2.7 (see Fig. QJ. 

One observes that the short range correlation length 
Ashort is 0.6 for peripheral collisions and decreases to- 
wards Aghort = 0.4 in central collisions. This decrease 
of the short rang e correlation length is consistent with 
Ref. 1^3, El E3| and was speculated to be a signal for 
the formation of a quark gluon plasma. The fraction of 
correlated particles decreases from almost 100% in pe- 
ripheral collisions to 88% for central reactions. 



IV. CENTRALITY DEPENDANCE OF a 2 c 

In Fig. we show the experimental results from the 
UA5 and the PHOBOS experiment on the fluctuations 
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FIG. 2: a 2 c as a function of Ar] for 0% - 20% central Au+Au 
at rj = 2. Black squares are PHOBOS data p i l2ll I27I |28 | . 
The line is the model calculation using equation JJJ. 
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FIG. 3: oq versus rapidity r\ for 40% — 60% peripheral Au+Au 
at Arj = 0.5. Black squares are PHOBOS data |2ll2ll27ll2^|. 
The line is the model calculation using equation 0. 
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FIG. 4: Oq versus Ar? for 40% - 60% peripheral Au+Au at 
J7 = 2. Black squares are PHOBOS data [H miHEl ■ The 
line is the model calculation using equation JSJ. 
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FIG. 5: Cluster multiplicity k versus the number of partici- 
pants. The last point represents 0%-3% Au+Au is computed 
by extrapolating our model parameters in Table [I] 



of the forward-backward correlation for pp (Au+Au) 
at different centralities. 

From Fig. El we see that Oq starts at 2.2 for pp and 
then increases to 2.7 for 40% — 60% peripheral and then 
goes down to 2.2 again for 0% — 20% central Au+Au. 
This behavior has important consequences that will help 
to understand the particles productions mechanism in 
nucleus-nucleus collisions compared to nucleon-nucleon 
collisions as we will discuss next. 

The increase in the dynamic fluctuations for 40%— 60% 
peripheral collisions compared to pp at the same energy 
might be attributed to the difference in the dynamics of 
the collisions between pp and nucleus- nucleus collisions. 
First let us consider the string picture for particles pro- 
duction in pp collisions. During such collisions a few 
strings between the quarks in both protons form and are 
subsequently fragmented into hadrons. If the same pic- 
ture holds for nucleus- nucleus collisions, then Oq in both 
pp and Au+Au will be similar and Oq has a flat centrality 
dependence. However, experimental data indicates that 
the cluster multiplicity k is increased in 40% — 60% pe- 
ripheral Au+Au collisions compared to pp interactions. 
An explanation of this enhancement might be found in 
the context of string fusion models |3^|. Here, the fu- 
sion process might increase the probability of high mass 
clusters compared to pp collisions, resulting in a higher 
multiplicity per cluster. 

If the formation of heavier clusters is the explanation 
for the increase of k in Au+Au compared to pp, then 
one expects an even further increase of k towards central 
Au+Au collisions. This enhancement will be observable 
if interactions with the surrounding medium are negligi- 
ble compared to the increase in k due to string fusion. 
Fig. 5 shows that this is not the case, because the clus- 



ter multiplicity in 0% — 20% central Au+Au collisions de- 
creases to k = 2.2 . In addition, the number of correlated 
particles / also decreases towards more central Au+Au 
reactions (cf. Table HJ) . This indicates that in central 
collisions the interactions with the hot and dense QCD 
matter results in smaller clusters (reducing the value of k) 
as compared to peripheral collisions. With even tighter 
centrality cuts one might be able to see a complete melt- 
ing of the clusters leading to independent uncorrelated 
particle emission patterns (i.e. = 1 and k = 1). 

To test this idea, we linearly extrapolate those param- 
eters in Table □ to 0% - 3% Au+Au collisions. The last 
point in Fig. [5] corresponds to this centrality cuts. We 
find that k ss 1.9, A s hort ~ 0.3. This means that clus- 
ters formed at the most central collisions are partially 
melted compared to mid-peripheral collisions. Also the 
short range correlation length A s hort becomes smaller. It 
is worth to mention that the decrease of the rapidity cor- 
relation length seems to be a general feature of nucleus- 
nucleus collisions at high energies since the width of the 
balance function expected to decrease with centrality as 
a result of late hadronization as mentioned in Ref. [34 . 
This decrease is observed in the data reported in p35ll36f . 

To distinguish between hadronic and QGP effects, Oq 
should be measured at different energies and ion sizes 
and also with different transverse momentum cuts to in- 
clude/exclude the effects of mini-jets, since jets hadronize 
into a shower of hadrons which might mimic a clustering 
effect (within the jet cone). Thus, limiting the measure- 
ment of <7p to the soft region (i.e. pt < 2 GeV) allows to 
minimize jet artifacts. 

In this case conclusive statement can be given by com- 
paring peripheral data to the central one to extract / and 
k for different initial conditions and to study how both / 
and k change as a function of the initial partonic volume, 
life time and energy density. 



V. SUMMARY AND CONCLUSION 

In this letter we present a model to analyze forward- 
backward dynamical fluctuations as recently measured 
by the PHOBOS collaboration in Au+Au reactions at 
•y/SNN = 200 GeV. We find that the effective cluster 
multiplicity in central collisions is lower than for periph- 
eral collisions. Also the short range rapidity correlation 
length decreases towards more central interactions. 

This might indicate that the clusters are "melted" in 
central collisions, resulting in a smaller fraction of corre- 
lated particles (decrease of /). One can test the degree 
of melting by measuring k as a function of centrality, and 
especially for very central Au+Au collisions. If one ob- 
serves that k — ► I, all clusters are melted and the final 
particles are independently emitted. 
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